program q2
	real	::	c, xmax, t
	integer	::	i, n
	real	:: x(100), y(100)

	c = 0.002
	xmax = sqrt(1.0/c)
	i = 1
	print*,"c=",c
	do t = 0.1, xmax, xmax/100.0
		print*, t, (1.0/(c**(2.0/3.0))) * t**(2.0/3.0) - t**2
	enddo
end program
